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We provide a short historic of the early development of kinetic theory in plasma physics and syn- 
thesize the basic kinetic equations describing the evolution of systems with long-range interactions 
derived in Paper I. We describe the evolution of the system as a whole and the relaxation of a test 
particle in a bath at equilibrium or out-of-equilibrium. We write these equations for an arbitrary 
long-range potential of interaction in a space of arbitrary dimension d. We discuss the scaling of 
the relaxation time with the number of particles for non-singular potentials. For always spatially 
homogeneous systems, the relaxation time of the system as a whole scales like in d > 1 and like 
A*'^ (presumably) in d = 1. For always spatially inhomogeneous systems, the relaxation time of the 
system as a whole scales like A'^ in any dimension of space. For one dimensional systems under- 
going a dynamical phase transition from a homogeneous to an inhomogeneous phase, we expect a 
relaxation time of the form A'''' with 1 < 5 < 2 intermediate between the two previous cases. The 
relaxation time of a test particle in a bath always scales like A'". We also briefly discuss the kinetic 
theory of systems with long-range interactions submitted to an external stochastic potential. This 
paper gathers basic equations that are applied to specific systems in Paper III. 

I. INTRODUCTION 

The kinetic theory of systems with long-range interactions is currently a topic of active research [H-Q. Initially, 
kinetic theories were developed for 3D plasmas and 3D stellar systems. However, they are now considered at a more 
general level for different types of potentials of interaction in different dimensions of space. It is therefore desirable 
to develop a general formalism valid for arbitrary potentials with long-range interactions and, from that formalism, 
consider different applications. This is the purpose of this series of papers. In Paper I Q, using the Klimontovich 
approach, we have derived the basic kinetic equations governing the evolution of spatially homogeneous systems with 
long-range interactions (these equations may be derived equivalcntly from the BBGKY hierarchy). In Paper HI 
we apply them to specific systems. In particular, we consider systems with power-law potentials (including plasmas 
and stellar systems in one, two, and three dimensions) and the Hamiltonian Mean Field (HMF) model with attractive 
or repulsive interactions. In order to facilitate the applications, we summarize in this paper the basic equations and 
provide additional results that were not given previously. Although these equations have been established for a long 
time in plasma physics for the 3D Coulombian interaction, no synthesis of these equations and of their generalization 
to other potentials exists yet. We think therefore that the synthesis given in this paper may be useful to many readers 
working in the field. We also provide a short historic of the early development of kinetic theory in plasma physics 
because it is in that context that the basic kinetic equations (Landau, Vlasov, Lenard-Balescu) were first derived. 
This historic is certainly incomplete but it contains references to important contributions that do not seem to be 
well-known. It may therefore help reconstructing in detail the body of knowledge of this topic. 

The first kinetic theory was developed by Boltzmann (1872) [9| for a dilute neutral gas in which the atoms or the 
molecules have short-range interactions. He obtained a kinetic equation by using simple semi-intuitive arguments. 
If we apply it to systems of particles interacting with power-law potentials decreasing at large distances as r~" we 
find that collisions with large impact parameters produce a divergence when a < 1 in d = 3. The critical exponent 
ttc = 1 precisely corresponds to the Coulombian interaction. Therefore, the kinetic theory of neutral gases developed 
by Boltzmann is not well-suited to handle inverse-square forces between ionized particles in a plasma. In the classical 
theory of neutral gases [loj , the assumption is made that only the relatively close encounters are important and that 
the forces between particles at greater distances have no effect. This is precisely the inverse situation that is met in 
a plasma. Distant encounters, producing small deflections, are more important than close encounters. It is therefore 
necessary to develop a specific kinetic theory for plasmas. 

The first kinetic equation of plasma physics was derived by Landau (1936) [Tl| who considered the effect of binary 
collisions between electrons and ions in a spatially homogeneous Coulombian plasma. He started from the Boltzmann 
equation and made a weak defiection approximation. Indeed, for the Coulomb potential, slowly decreasing with the 
distance like r~^, weak collisions are the most frequent ones (the probability of large scattering is small). Each 
encounter induces a small change in the velocity of a particle but the cumulated effect of these encounters leads to 
a macroscopic process of diffusion and friction that drives the system towards statistical equilibrium. It is therefore 
possible to expand the Boltzmann equation in terms of a small deflection parameter A <C 1 and make a linear 
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trajectory approximation. This treatment yields a logarithmic divergence at small and large impact parameters but 
the equation can still be used successfully if appropriate cut-offs are introduced. The divergence at small scales 
arises because the Landau theory, based on a weak coupling approximation, docs not take proper account of the 
processes occurring in close binary collisions, where the trajectories of the particles are no more rectilinear. A natural 
lower cut-off Xl, which is called the Landau length, corresponds to the impact parameter leading to a deflection at 
90°. The divergence at large scales arises because the Landau theory does not take collective effects into account. 
Phenomcnologically, the Debye length^ provides a natural upper cut-off. Indeed, beyond the Dcbye length Ad many 
particles interact and the electric field of a particle is screened by its neighbors. Therefore, the particles do not interact 
if their separation exceeds the Debye length. In this sense, the interaction is effectively short-ranged. 

Since the Coulombian interaction between charges shares many analogies with the Newtonian interaction be- 
tween masses, it is useful to discuss in parallel the kinetic theory of self-gravitating systems. Chandrasekhar 
(1942, 1943a, 1943b) [T^ - [l6j developed a kinetic theory of stellar systems in order to determine the timescale of colli- 
sional relaxation and the rate of escape of stars from globular clusters^. To simplify the kinetic theory, he considered 
an infinite and homogeneous system and studied the relaxation of a test particle in a thermal bath. He started from 
the Fokker-Planck equation and determined the diffusion coefficient and the friction force (second and first moments 
of the velocity increments) by considering the mean effect of a succession of two-body encounters. This approach 
was based on Jeans (1929) [1^ demonstration that the cumulative effect of the weak deflections resulting from the 
relatively distant encounters is more important than the effect of occasional large deflections produced by relatively 
close encounters. Since his approach can take large deflections into account, there is no divergence at small impact 
parameters and the gravitational analogue of the Landau length appears naturally in the treatment of Chandrasekhar. 
However, his approach leads to a logarithmic divergence at large scales that is more difficult to remove in stellar dy- 
namics than in plasma physics because of the absence of Debye shielding for the gravitational force. In a series of 
papers, Chandrasekhar and von Neumann (1942) [2^ developed a completely stochastic formalism of gravitational 
fluctuations and showed that the fluctuations of the gravitational force are given by the Holtzmark distribution (a par- 
ticular Levy law) in which the nearest neighbor plays a prevalent role. From these results, they argued (see also Jeans 
[l9| and Spitzer [2lj) that the logarithmic divergence has to be cut-off at the interparticle distance I ^ n~^^^ where n 
is the numerical density of stars. However, since the interparticle distance is smaller than the Debye length, the same 
arguments should also apply in plasma physics, which is not the case. Therefore, the conclusions of Chandrasekhar 
and von Neumann are usually taken with circumspection. In particular, Cohen et al. (1950) [2^ argued that the 
logarithmic divergence should be cut-off at the Jeans length since the Jeans length is the gravitational analogue of the 
Debye length. Indeed, while in neutral plasmas the effective interaction distance is limited by the Debye length, in a 
self-gravitating system, the distance between interacting particles is only limited by the system's size. Chandrasekhar 
(1949) [2^ also developed a Brownian theory of stellar systems and showed that, from a qualitative point of view, 
the results of kinetic theory can be understood very simply in that framework as a competition between diffusion 
and friction'^. In particular, he showed that a dynamical friction is necessary to maintain the Maxwell-Boltzmann 
distribution of statistical equilibrium and that the coefficients of friction and diffusion are related to each other by 
an Einstein relation which is a manifestation of the fluctuation-dissipation theorem. This relation is conflrmed by 
his more precise kinetic theory based on two-body encounters. The Fokker-Planck approach of Chandrasekhar was 
further developed by Rosenbluth et al. (1957) [l^l who proposed a simplified derivation of the coefficients of diffusion 
and friction valid for an arbitrary distribution function. It is important to emphasize, however, that Chandrasekhar 



^ The Debye length was first introduced in connection with screening processes in highly ionized electrolytes by Debye and Huckcl (1923) 
I12I Il3l| . Each particle of the plasma tends to attract to it particles of opposite charge and to repel particles of like charge, thereby 
creating a kind of "cloud" of opposite charge. This cloud screens the electric field produced by the particle on a distance of the order 
of the Debye length. 

^ Early estimates of the relaxation time of stellar systems were made by Sehwarzschild (1924) fl7| . Rosseland (1928) [T3 |. Jeans (1929) 
Tsl . Smart (1938) |20(I . and Spitzer (1940) |21| . On the other hand, the evaporation time was first estimated by Ambartsumian (1938) 
H and Spitzer (1940) [H. 

^ The stochastic evolution of a star is due primarily to many small deflections produced by relatively distant encounters. The problem 
of treating particles undergoing numerous weak deflections was originally encountered in relation to the Brownian motion of large 
molecules which are thermally agitated by the smaller field molecules. The stochastic character of the many small impulses which act on 
a suspended particle is usually described by the Fokker-Planck equation. Chandrasekhar [2^1 (1943) noted the analogy between stellar 
dynamics and Brownian theory and employed a Fokker-Planck equation for the velocity distribution function. Indeed, he argued that a 
given star undergoes many small-angle (large impact parameter) collisions in a time small compared with that in which its position or 
velocity changes appreciably. Since the Newtonian potential makes the cumulative effect of these small momentum transfers dominant, 
the stochastic methods of the Fokker-Planck equation should be more appropriate for a stellar system than the Boltzmann approach. 
However, the coefficients of diffusion and friction are finally calculated with the binary-collision picture. This is equivalent to expanding 
the Boltzmann equation for weak defiexions (as in the Landau theory) since the Fokker-Planck equation can precisely be obtained from 
the Boltzmann equation in that limit. 
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did not derive the kinetic equation for the evolution of the system as a whole. Indeed, he considered the Brown- 
ian motion of a test star in a fixed distribution of field stars (bath) and derived the corresponding Fokker-Planck 
equation"*. King (1960) [l^ noted that if we were to describe the dynamical evolution of the cluster as a whole, the 
distribution of the field particles should evolve in time in a self-consistent manner so that the kinetic equation must 
be an integrodifferential equation. The kinetic equation obtained by King turns out to be equivalent to the Landau 
equation, although written in a different form^. There is, however, an important difference between stellar dynamics 
and plasma physics. Neutral plasmas are spatially homogeneous due to clectroncutrality and Dcbye shielding. By 
contrast, stellar systems arc spatially inhomogencous. The above-mentioned kinetic theories developed for an infinite 
homogeneous system can be applied to an inhomogencous system only if we make a local approximation. In that case, 
the collision term is calculated as if the system were spatially homogeneous or as if the collisions could be treated as 
local. Several astrophysicists have then tried to incorporate effects of spatial inhomogeneity in the kinetic theory of 
stellar systems (see further discussion in jsoj). 

Returning to plasma physics, the early works on kinetic theory were based on the model of binary collisions 
treated either by the Boltzmann or, more adequately, by the Fokker-Planck equation. However, these theories lead to 
difficulties because the collision term diverges at large impact parameters. This is due to the fact that the Coulomb 
potential docs not decrease rapidly enough at large distances, which means physically that collective interactions, 
involving many particles simultaneously, are important in determining the behavior of the system. In plasma physics, 
this divergence has been suppressed by introducing rather arbitrarily a cutoff in the integration at the Dcbye length 
which roughly measures the effective screening of the interaction by the polarization of the medium®. However, this 
ad hoc procedure is not satisfactory and a proper treatment of collective effects must be developed. Indeed, the 
encounters are clearly not binary since there are many particles inside the screening sphere. 

A first treatment of collective effects in plasma physics was introduced by Vlasov (1938,1945) [s^, H^]. He treated 
collective effects in a self-consistent manner but neglected "collisions" between charges. Therefore, each charge is 
described as moving under the infiuence of an average field due to all the others. The evolution of the distribution 
function is then governed by the coUisionless Boltzmann equation coupled to the Poisson equation. This system of 
equations had been previously introduced by Jeans (1915) [38| for stellar systems under the same approximations (see 
Henon (1982) ^] for interesting historical details). However, the Vlasov equation cannot be used for a description 
of transport processes. Because of the neglect of collisions, it contains no mechanism of relaxation toward thermal 
equilibrium . 

Pines and Bohm (1952) [42| developed for the first time a description of electron interactions taking collective 
effects and binary collisions into account. They showed that for phenomena involving distances greater than the 
Debye length, the system behaves collectively while for distances shorter than the Debye length, it may be treated 
as a collection of approximately free individual particles, whose interaction may be described in terms of two-body 
collisions. This justifies therefore the Landau procedure of cutting-off the integration at the Debye length. They 
studied screening properties of out-of-equilibrium plasmas. For particles at rest or moving sufficiently slowly, they 
showed that the field is screened as in the Debye theory. However, if particles are moving with a speed greater than 
the mean thermal speed, the field propagates away from the source as a wake. As a result, an excess of particles 
of like-charge accumulates ahead of the moving particle and conversely an excess of particles of opposite-charge 
accumulates behind it. Under such circumstances, polarization occurs and the potential field becomes asymmetrical. 
They concluded that the interaction between electrons in a plasma can be described in terms of short-range collisions 
between "effective free particles" corresponding to the electrons plus their associated cloud. 

Bogoliubov (1946) [43] was the first to improve the Landau kinetic theory by taking collective effects into account. 



* Indeed, Chandrasekhar Il5l . ITsIl models the evolution of globular clusters by the Kramers equation which has a fixed temperature 
(canonical description) while a more relevant kinetic equation would be the Landau equation which conserves the energy (microcanonical 
description). 

^ App arently, the Landau equation was not known by astrophysicists at that time who rather followed the treatment of Chandrasekhar 
Il4l4l6ll and Rosenbluth et al. |27ll . To our knowledge, the first explicit reference to the Landau equation in the astrophysical literature 
appeared in the paper of Kandrup (1981) [29II . There is also a strange comment related to the work of Landau in the paper of Cohen 
et al. (1950) [2^ that sug gests that his work was not clearly understood. 

^ Before Landau (1936) [ll| . several authors like Gruner (1911) [3l| . Chapman (1922) [3^ . and Persico (1926) [ssj attempted to develop a 
kinetic theory of ionised gases to model the interior of giant stars. Gruner [3l| applied the kinetic theory of gases to power-law potentials 
r ' and noted that a divergence occurs when 7 = 1. Chapman |32|| proposed to truncate the diverging integral at the interparticles 
distance (the Debye shielding was not yet discovered) . Persico [33^ , following the works of Rosseland |34{ and Eddington |35| | , understood 
that the divergence could be solved by taking collective effects (Debye shielding) into account. He replaced the Coulombian potential 
of interaction by the Debye potential in the expression of the diffusion coefficient coming from the theory of gases. 

^ There is, however, an irreversible phenomenon of Landau damping (1946) [40| which is a process of uniformization in configuration space. 
In stellar dynamics, the Vlasov- Poisson system may also undergo a process of violent coUisionless relaxation that has been discussed by 
Lynden-Bell (1967) [4l| in terms of statistical mechanics. 
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Starting from the Liouville equation and using what is now called the "BBGKY hierarchy" and the "Bogoliubov 
ansatz" , he derived a kinetic equation in which the collision term is expressed in terms of a two-body correlation 
function that is solution of an integral equation (this equation corresponds to the first term of an expansion in powers 
of the coupling parameter 5 = 1/ (nX'j^) so it is valid in a "weak coupling" approximation). In this equation, the Debye 
shielding is automatically taken into account so that the integral over the space coordinates is convergent at large 
distances contrary to the treatment of Landau who had to introduce a cutoff artificially at the Debye length in order 
to make a logarithmically divergent integral finite. In addition, while the Boltzmann and Fokkcr-Planck equations 
used by Landau and Chandrasekhar have a phenomcnological character, the method of Bogoliubov, that starts from 
the Liouville equation, is systematic and rigorous. However, the final equation for the evolution of the distribution 
function given by Bogoliubov is not explicit. 

Lenard (1960) [31 solved the Bogoliubov integral equation exactly by using a Fourier transformation and obtained 
an explicit kinetic equation where the time derivative of the distribution function is expressed explicitly in terms 
of the distribution function itself. The same kinetic equation, describing the collective effects in a weak coupling 
approximation, was obtained independently by Balescu (1960a, 1960b) [1^,14^1 from the general theory of irreversible 
processes in gases developed by Prigogine and Balescu (1957) [l^. In this formalism, one starts from the Liouville 
equation, perform a Fourier transformation and solve the obtained equations by an iteration procedure. The asymp- 
totic evaluation of the orders of magnitude of the various contributions to the Fourier components is accomplished by 
means of a diagram technique. This leads to an integral equation equivalent to the Bogoliubov integral equation that 
Balescu solved by using essentially the same method as Lenard. This is valid in a weak coupling approximation at the 
order g = 1/A where A = nX'jj is the number of charges in the Debye sphere. The Lenard-Balescu kinetic equation 
describes the interaction of "quasi particles" , which are electrons or ions "dressed" by their polarization clouds (this 
picture in terms of "effective free particles" -electron plus its associated cloud- was previously suggested by Pines and 
Bohm). These clouds are not a permanent feature, as in equilibrium theory, but they have a nonequilibrium changing 
shape distorted by the motion of the particles. The Lenard-Balescu equation exhibits a new type of nonlinearity 
which is directly related to the collective nature of the interactions. Besides the product of two distribution functions 
/(l)/(2) characteristic of any two-body collision term, it contains the distribution function in the denominator. In 
fact, the Lenard-Balescu equation differs from the original Landau equation by the appearance of the dielectric func- 
tion |e(k, k • v)| in the denominator of the potential of interaction u(k). Physically, this means that the particles are 
"dressed" by their polarization cloud. The Landau equation is recovered when |e(k, k • v)| is replaced by unity, i.e. 
when collective effects are neglected. However, with this additional term, it is found that the logarithmic divergence 
that occurs at large scales in the Landau equation is removed and that the Debye length appears naturally. In the 
dominant approximation, the Lenard-Balescu equation reduces to the Landau equation with a large-scale cut-off at 
the Debye length (which is now rigorously justified). 

The Lenard-Balescu equation describes the interaction of the electrons which are "dressed" by a polarization 
could^. This effect is well-known in equilibrium as the Debye cloud. However, in the situation considered by Lenard 
and Balescu, the polarization cloud differs sensibly from the latter because the system is out-of-equilibrium and 
the polarization cloud is deformed by the motion of the particles. For small velocities the deformation is merely a 
fiattening in the direction of the motion but for velocities larger than the thermal velocity the cloud degenerates 
into a wake of plasma oscillations. This precisely accounts for what Pines and Bohm have shown in a more heuristic 
manner. For small velocities, or in the dominant approximation, the effect of screening in the Lenard-Balescu equation 
is equivalent to replacing the bare Coulomb potential by the Debye potential (in agreement with Pines and Bohm). 
Interestingly, Liboff (1959) [i^ had previously developed a kinetic theory of plasmas with a Debye potential in order 
to avoid divergences at large scales^. He showed that, in the dominant approximation, the Debye potential gives the 
same results as the "Coulomb plus Debye cutoff" . 

The Lenard-Balescu equation was also derived independently by Guernsey (1960) so it is sometimes called the 
Lenard-Balescu-Guernsey equation^". The works of Lenard, Balescu and Guernsey have removed the problematic 
divergence at large scales that appears in the Landau equation. However, their approach based on a weak coupling 
approximation yields a divergence at small scales because it does not take strong collisions into account. Coulombian 
plasmas therefore seem to be described by two complementary kinetic theories, depending on the value of the impact 
parameter A: (i) the binary encounters (or impact) theory based on the Boltzmann or Fokker-Planck equations which 
is appropriate to describe strong collisions when A ~ Al (it does not present any divergence at small scales but 
yields a logarithmic divergence at large scales) and (ii) the wave theory based on the Lenard-Balescu equation that is 



* This notion of "dressed test particles" is further developed by Rostoker (1964) [48l |. 

^ Liboff computes the diffusion and friction coefficients from the two-body encounters theory taking strong collisions into account, so he 

has no divergence at small scales neither. 
" The kinetic theory of Guernsey, and other kinetic theories, are presented in the book of Wu (1966) [sill . 
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appropriate to describe collective effects when X ^ Xd (it does not present any divergence at large scales but yields a 
logarithmic divergence at small scales). In the intermediate case <^ A ^ A^, these two kinetic theories reduce to 
the Landau theory (it presents divergences both at small and large scales). Baldwin (1962) [s^], Frieman and Book 
(1963) [l^l, Weinstock (1964) [s^l and Guernsey (1964) proposed a kinetic equation that takes strong collisions and 
collective effects into account and exhibits no divergence. Their treatment essentially amounts to writing the collision 
term as a sum of two terms: a Lenard-Balescu term that takes collective effects into account and a Boltzmann term 
that takes strong collisions into account. This leads to a general equation combining the Boltzmann and the Lenard- 
Balescu equation. Another unification scheme has been proposed by Kihara and Aono (1963) |5Q] by introducing an 
artificial splitting of the integration over the impact parameter in the impact theory and over the wavenumber in the 
wave theory. Finally, Gould and DeWitt [s^l (1967) obtained a convergent kinetic equation by using a direct analogy 
with the equilibrium theory of the electron gas. They also carefully compared their results with those of previous 
works. 

Prior to Lenard and Balescu, several authors like Temko (1957) [Hi, Kadomtsev (1958) Tchen (1959) 
Ichikawa (1960) [6l| and Willis (1962) had attempted to solve the Bogoliubov integral equation. However, they 
solved it in an approximate manner which amounts essentially to making the Debye-Hiickel approximation, i.e. to 
using the equilibrium form of the correlation function. Although their treatment is approximate, it shows that the 
correlations that exist in equilibrium are sufficient to provide both short and long distance cutoffs. Furthermore, 
the treatment of Tchen goes beyond the Debye-Hiickel approximation and captures the ellipsoidal anisotropy of the 
shielding cloud first evidenced by Pines and Bohm for a moving particle. These authors obtained explicit kinetic 
equations that have the form of the Landau equation in which the Coulombian potential of interaction is replaced by 
the Debye-Hiickel potential (or more general potentials). This justifies the procedure of Liboff who introduced the 
Debye-Hiickel potential phcnomcnologically. In the dominant approximation, the velocity dependence of the dielectric 
function arising in the Lenard-Balescu equation may be neglected, so these heuristic treatments arc sufficient for 3D 
plasmas. 

In parallel, several authors like Gasiorowicz et al. (1956) l63l . Temko (1957) [H^, Tchen fl959') [60j. Rostoker 
and Rosenbluth (1960) [H, Thompson and Hubbard (1960) jfill, and Hubbard (1960a,1960b) [H, [ejdeveloped a 
kinetic theory of homogeneous plasmas for a test particle in a thermal bath, taking collective effects into account. 
They showed that the large-scale divergence that appears in the Landau diffusion coefficient and friction force is 
actually regularized at the Debye length. In particular, Hubbard [g^, [63] considered the relaxation of a test particle 
in a thermal bath (one particle out-of-equilibrium moving in a gas that is in equilibrium), or in a bath that is not 
necessarily in thermal equilibrium, and derived the diffusion and friction coefficients entering in the Fokker-Planck 
equation. He showed that the diffusion is due to the fluctuations of the electric field while the friction is due to a 
polarization process^^. The field produced by the test particle polarizes the plasma. The charge displacement in the 
plasma representing this polarization is highly correlated with the position of the test particle and gives rise to an 
electric field which acts on the test particle and produces a friction. His calculations take into account the effect of 
correlations between distant electrons or ions. Consequently, the Debye screening is automatically included in his 
kinetic theory and it is not necessary (as it had been in the earlier calculations) to introduce a cut-off procedure to 
suppress the effects of distant encounters between particles. His calculations also include a proper treatment of close 
binary encounters like in the standard Boltzmann theory described by Chapman and Cowling [lo| . As a result, he 
obtains Fokkcr-Planck coefficients which include both collective effects and the contribution of close binary encounters. 
In his theory, no ad hoc cut-off procedures of any kind are needed. In the dominant approximation, he recovers the 
results of Landau but now without any divergence. The Landau and Debye lengths appear naturally in his treatment. 
It is interesting to note, for historical reasons, that the results obtained by Hubbard are very closely related to those 
obtained independently by Lenard and Balescu at the same time. Indeed, if we substitute the diffusion and friction 
coefficients found by Hubbard in the Fokker-Planck equation and perform minor transformations (a substitution that 
Hubbard did not explicitly make), one obtains the Lenard-Balescu equation! 

In the following years, alternative derivations of the Lenard-Balescu equation have been given. Some derivations 
are based on the BBGKY hierarchy but use a procedure different from that of Lenard and Balescu. For example, 
Dupree [11] and Ichimaru [g^ solve the second equation of the BBGKY hierarchy in terms of operators (propagators) 
which are related to the solution of the linearized Vlasov equation. On the other hand, Nicholson [70[ solves the 
second equation of the BBGKY hierarchy by using Fourier-Laplace transforms (while Lenard and Balescu use only 
Fourier transforms). Other derivations due to Dupree [7l[, Fried Klimontovich [t^], and Lifshitz and Pitaevskii 



This interpretation is particularly clear in the work of Gasiorowicz et al. [63| . As a particle moves through the distribution of particles, 
the latter becomes polarized because the test particle attracts the field particles, so that they tend to concentrate behind it, and thus 
the force from particles behind is larger than the force due to those in front, with the result that the test particle is slowed down. 
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[TJ] start from the Klimontovich equation and use a quasilinear approximation. All these methods are equivalent 
but the calculations are presented in a different manner. It is usually acknowledged that the approach based on the 
Klimontovich equation is technically simpler. 

Recently, it was realized by the community working on systems with long-range interactions that these kinetic 
equations have a wider scope than just plasma physics and stellar dynamics. Actually, they are valid for arbitrary 
potentials of long-range interactions, attractive or repulsive, in a space of dimension d. Although the form of the 
Lenard-Balescu equation is unchanged (we just need to replace the Fourier transform of the Coulombian potential by 
the Fourier transform of the more general potential w(fc)), the results (relaxation time, diffusion coefficient, polarization 
cloud...) sensitively depend on the form of the potential of interaction and on the dimension of space. It is therefore 
useful to expose general kinetic equations valid for arbitrary potentials with long-range interactions, then consider 
different applications. 



II. THE THERMODYNAMIC LIMIT AND THE WEAK COUPLING APPROXIMATION 

We consider a system of material particles in interaction in a d-dimensional space. Their dynamics is described by 
the Hamilton equations 

dri dH dvi dH 

dt dvi ' dt dvi ' 

with the Hamiltonian 

2—1 i<j 

where u(|r — r'|) is the potential of interaction and m the individual mass of the particles. The first term K in the 
Hamiltonian is the kinetic energy and the second term U is the potential energy. We assume that the potential of 
interaction decays at large distances as r~'^ with j < d. In that case, the potential is said to be long-ranged. For 
such potentials, the energy J^°° u{r)r'^~^ dr diverges at large scales implying that all the particles interact with each 
other and that the system displays a collective behavior. 

We introduce a characteristic length R and a characteristic velocity Vm- We define the dynamical time hy to = 
R/Vm- For systems with long-range interactions, the potential energy scales as C/ ^ N^m^u{R) while K ^ Nmv^. The 
kinetic energy and the potential energies are comparable if Nmv^ ^ N^m^u{R) which implies that ^ Nmu{R). 
As a result, the energy scales as -E ~ Nmv^ ^ N'^m?u{R) and the kinetic temperature, defined by kgT ~ mw^, 
scales as ksT ^ Nm?u{R) ^ E/N . Inversely, these relations may be used to define R and Vm as a function of the 
energy E (conserved quantity in the microcanonical ensemble) or as a function of the temperature T (fixed quantity 
in the canonical ensemble). The proper thermodynamic limit of systems with long-range interactions corresponds to 
N — )■ -t-oo in such a way that the normalized energy e = E/N'^m'^u{R) and the normalized temperature r] = f3Nm,'^u{R) 
are of order unity. We introduce the coupling parameter g = Epot/Ekm = m?u{R) jkBT where Epot ~ m^u{R) is the 
potential energy of two particles separated by the distance R and Ekin ^ ksT is the thermal energy. According to 
the previous estimates, we have g ^ 1/N. Therefore, when N — > +oo, we can consider a weak coupling approximation 
since g 1. 

It is convenient to rescale the distances by R, the velocities by Vm, the time hy tjj, and the mass by m. This is 
equivalent to taking R = Vm = tz)=m = lm. the original equations. In order to satisfy the condition Nmv'^ ^ 
N^m^u{R) the potential of interaction must scale like u{R) ~ It is therefore convenient to write u{r) = jju{r) 

with u{R) ^ 1. This is the Kac prescription for long-range interactions [zHl- With this normalization, we have E ^ N, 
S N and T ~ 1 in the limit N — >■ -|-oo. The energy and the entropy are extensive but they remain fundamentally 
non-additive The temperature is intensive. This normalization is very convenient since the length, velocity, time 
and mass scales are of order unity. Furthermore, since the coupling constant u scales as 1/-/V, this immediately shows 
that a regime of weak coupling holds when N ^ 1. 

From the Liouville equation, we can derive the BBGKY hierarchy of equations for the reduced distribution functions. 
This hierarchy of equations can also be expressed in terms of the correlation functions (an equivalent hierarchy of 
equations may be obtained from the Klimontovich equation) . The previous rescaling shows that only the dimensionless 
(coupling) parameter g ~ 1/N arises in the equations. For N — > +oo, we can consider an expansion of the equations 
of the BBGKY hierarchy in powers of this small parameter g = <C 1. This corresponds to the weak coupling 
approximation. It is usually argued that the correlation functions scale as P'^ ^ 1/N"~^. In particular, the two-body 
correlation function scales as ~ 1/^ ^^'^ the three-body correlation function scales as ~ l/N"^. This gives a 
systematic procedure to expand the BBGKY hierarchy in powers of the coupling parameter g = 1/N <^ 1. 
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For self-gravitating systems in which the constituents interact with a potential u = ~G/r, we introduce the Jeans 
wavenumber kj = (47rG/3mp)^/^ and the gravitational pulsation ojg = {4:nGpY^'^ , where p = nm is the mass density 
and /3 = l/{kBT) the inverse kinetic temperature. We may use the Jeans length Aj ~ 27r/fcj as a relevant lengthscale 
and the dynamical time to = 2'k/ujq as a relevant timescale. From the virial theorem, the Jeans length gives an 
estimate of the system's size R and the dynamical time may be written as to ^ R/vm- We note that these scalings 
may be obtained from the general arguments given above. We then find that the only dimensionless parameter 
in the problem is the coupling parameter g = \/K where A = riAj ~ N gives the typical number of particles in 
the Jeans sphere (i.e. in the system). Alternatively, if we define units of size, time, velocity and mass such that 
A J = to = v„i = m = 1, we must take G ^ 1/N for consistency. The weak coupling approximation corresponds to 
A iV ^ +00. 

For Coulombian plasmas in which the constituents interact with a potential u = (e^/m^)/r, we introduce the 
Debye wavenumber kjj = (47re^/?p/m)^/^ and the plasma pulsation wp = (47re^p/m^)^/^, where p = nm is the mass 
density. We may use the Debye length A_d = 2Tr/ kjj as a relevant lengthscale and the dynamical time to = 27r/wp 
as a relevant timescale. The Debye length gives an estimate of the effective range of interaction due to screening by 
opposite charges and the dynamical time may be written as ~ Xo/vm- We then find that the only dimensionless 
parameter in the problem is the coupling parameter g = 1 /A where A = nX^ gives the typical number of particles in 
the Debye sphere. Alternatively, if we define units of size, time, velocity and mass such that Xjj = tjj = = m = I), 
we must take ~ 1/A for consistency. The weak coupling approximation corresponds to A — )■ +00. 

We show that 1/A may indeed be interpreted as a coupling parameter. The coupling parameter F is defined as the 
ratio of the interaction strength at the mean interparticle distance Grn^v}^^ (resp. e^'n}/^) to the thermal energy fc^T. 
This leads to F = Gm^n^l'^ jkeT = XjijiX^jfl'^ = l/A^/s ^ 1/^2/3 (j.ggp_ y = e^n^/^/kBT = l/{nXl,fl'^ = l/A^/S). 
If we define the coupling parameter g as the ratio of the interaction strength at the Jeans (resp. Debye) length Gm} jXj 
(resp. e^/Ap))) to the thermal energy fc^T, we get g = 1/A. Therefore, the expansion of the BBGKY hierarchy in 
terms of the coupling parameter F or g is equivalent to an expansion in terms of the inverse of the number of particles 
in the Jeans sphere A = riA j ^ N (resp. the inverse of the number of particles in the Debye sphere A = nA^,). The 
weak coupling approximation is therefore justified when A ^ 1. 

The kinetic equations of systems with long-range interactions are well-known. They may be derived either from 
the BBGKY hierarchy or from the Klimontovich equation. For g = 1/N = 0, we get the Vlasov equation. This 
describes a "collisionless" evolution. In that case, the mean field approximation is exact. At the order g = 1/A^, we 
get the Lenard-Balescu equation. This is the first order correction to the Vlasov equation due to "collisions" (more 
properly correlations) between particles. The generalization of kinetic theory to spatially inhomogeneous systems, 
using angle-action variables, is more recent (30l l76l - l79| . The kinetic theory of 2D point vortices, which form a peculiar 
system with long-range interactions, is also relatively new (see [soj and references therein). In the following, we 
summarize the basic kinetic equations of spatially homogeneous systems with long-range interactions appropriately 
generalized to arbitrary potentials of interaction (see Paper I for the details of their derivation) and we discuss the 
results depending on the dimension of space. We distinguish the evolution of the system as a whole and the relaxation 
of a test particle in a bath that may be at equilibrium or out-of-equilibrium. We discuss the scaling of the relaxation 
time with the number of particles^^. In Paper III, we apply these kinetic equations to particular systems of physical 
interest. 



III. EVOLUTION OF THE SYSTEM AS A WHOLE: THE LENARD-BALESCU EQUATION 

A. Spatially homogeneous systems in d > 1 

Spatially homogeneous systems with long-range interactions are governed by the Lenard-Balescu equation 

I - "P"'"'!- / *<'^' heEf*''' ^ <^ - ^''1 ilk - ^) ^'^•"^<^' P> 

where u{k) is the Fourier transform of the potential of interaction and e(k, w) is the dielectric function 

r k- ^ 

e(k, w) = 1 + {2n f u{k) / ^ dv. (4) 

J w — k • V 



^ In this paper, we assume that the potential is non-singular at the origin. The case of pure power-law potentials that are divergent at 
the origin is considered in Paper III. 
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The distribution function /(v,t) is normahzed such that J f dv = p gives the mass density. This kinetic description 
assumes that the distribution function /(v,i) is Vlasov stable at any stage of the dynamics. The Lenard-Balescu 
equation is valid at the order 1/iV in the limit +00 (as explained in Sec. HIl this scaling may be seen on Eq. 

([5]) by using dimensionless variables such that |r| ~ 1, |v| ~ 1, to ^ 1, m ^ 1/-^, and working with the one-body 
distribution Pi = f /N ~ 1). This corresponds to the weak coupling approximation in which three-body collisions 
and strong collisions may be neglected. The evolution of the system is then due exclusively to weak collisions and 
collective effects (see Appendix A of Paper III). 

If we neglect collective effects, which amounts to replacing |e(k, k • v)| by 1, we get the Landau equation 



^ = ^(2^)'"^^ / d\^dv' hk,u{krS[k ■ (v - v')] - ^) /(v,i)/(v',i). 



(5) 



Inversely, the Lenard-Balescu equation may be obtained from the Landau equation by replacing the "bare" potential 
of interaction u{k) by a "dressed" potential of interaction 

, , u(k) 

Udressed{k, k • v) = — — j . (6) 

|e(k,k-v)| 

Physically, the particles are "dressed" by their polarization cloud. We note that the Lenard-Balescu equation takes 
into account dynamical screening since the velocity of the particles explicitly appears in the dressed potential of 
interaction. If we replace |e(k, k • v)| by |e(k, 0)| and approximate / by the Maxwellian distribution (p4|. we get the 
Dcbye-Hiickel potential (see Paper I): 

^DH (fc) = , , ,o,^L/L^«^. ■ (7) 



1 + {2T:Yu{k)l3mp' 



In that case, we obtain 



= TT(2TTfm-^ [ dkdVhk, ^(5[k • (v - v')] ( ^--^]/(v,t)/(v',i). (8) 



We shall call this equation the Debye-Hiickel approximation of the Lenard-Balescu equation. The Debye-Hiickel 
potential takes into account static screening (see Appendix \^ . Indeed, the limit w — > corresponds to the limit 
t — > +00 for the inverse Laplace transform of the effective potential ([6]). The static expression is reached when 
all the transient oscillations have been damped. Alternatively, the Debye-Hiickel potential (O may be seen as an 
approximation of Eq. ^ for small velocities. The Landau equation evaluated with the Debye-Hiickel potential will 
be an acceptable approximation of the Lenard-Balescu equation whenever the large- velocity population does not play 
a crucial role as compared to the bulk of the distribution. The validity of the Debye-Hiickel approximation of the 
Lenard-Balescu equation will be discussed in Paper HI. 

For the Landau equation ([5]) the integral over k can be performed explicitly (see, e.g., [T^I) and we get 

where w = v — v' is the relative velocity and Kd is a constant with value A'3 — Sn^m J^°° k^{i{k)^ dk in d — 3 and 

K2 = Stt^to /q^°° k^u{ky dk in d = 2. In the Debye-Hiickel approximation, Kd is replaced by Kj^^ defined in terms 
of UDH{k). The Landau equation may be expressed in terms of the Rosenbluth potentials as shown in Appendix [Bl 
The Lenard-Balescu equation ([3]) is valid at the order 1/iV so it describes the "collisional" evolution of the system 
on a timescale ~ Nto, where is the dynamical time. This kinetic equation conserves the mass M ~ J f dv and 

the energy E = J dv. It also monotonically increases the Boltzmann entropy S' = — / ^ In ^ dv. More precisely, 
for d > 1, we can show that S" > and 5" = if, and only if, / is the Maxwell-Boltzmann distribution (i/-theorem). 
The collisional evolution of the system is due to a condition of resonance between the orbits of the particles. The 
condition of resonance, encapsulated in the (5-function, corresponds to k • v' = k • v with v' 7^ v. For d > 1, the 
condition of resonance can always be satisfied. Therefore, the only steady state of the Lenard-Balescu equation is 
the Maxwell-Boltzmann distribution. This is the maximum entropy state at fixed mass and energy. Because of the 
H-theorem, the Lenard-Balescu equation relaxes towards the Maxwell-Boltzmann distribution for t — > +00. Since the 
collision term is valid at the order 0{1/N), the relaxation time scales as 

^whoie ^ j^^^^ (-^ ^ ^ homogeneous). (10) 
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B. Spatially homogeneous systems in d = 1 

For one-dimensional systems, the Lenard-Balescu kinetic equation ([3]) reduces to 

f - - (I; - 1:) /(... «/(./. « ^ 0, (11) 

where we have used the identity S{Xx) = jy^6{x). Therefore, the colhsion term vanishes at the order 1/7V because 
there is no resonance. The kinetic equation reduces to df /dt ^ so the distribution function docs not evolve at all on 
a timescale ~ Nt]j. To first order in l/N, the detailed balance between drag and diffusion is valid not only at thermal 
equilibrium but for any Vlasov stable distribution f(v). Therefore, the kinetic theory predicts no thermalization to a 
Maxwellian at first order in 1/iV. The Maxwellization is at least a second order effect in and consequently a very 
slow process. This result has been known for a long time in ID plasma physics [8l|, [s^] (in that context N represents 
the number of particles in the Debyc segment usually denoted A) and was rediscovered recently in the context of 
the HMF model (S^H^. This implies that, for onc-dimcnsional homogeneous systems, the relaxation time towards 
statistical equilibrium is larger than Ntjj. Therefore 

^whoie ^ j^^^^ ^ ^ homogeneous). (12) 

Since the relaxation process is due to more complex correlations than simply two-body collisions, we have to develop 
the kinetic theory at higher orders (taking three-body, four-body,... correlation functions into account) in order to 
obtain the relaxation time. If the collision term does not vanish at the next order of the expansion in powers of 1/iV, 
the relaxation time is of order 

^whoie _ ^2^^^ ^ ^ homogeneous). (13) 

This quadratic scaling, conjectured in [sij . was numerically observed in (85l.l86j for spatially homogeneous ID plasmas. 
In that case, the relaxation is caused by three-body correlations. 



C. Spatially inhomogeneous systems 

For spatially inhomogeneous systems, the Lenard-Balescu equation and the Landau equation must be written with 
angle-action variables [76l - l79l | . It is then found that spatial inhomogeneity allows for additional resonances [tJ . As a 
result, the relaxation time scales as 

^whoie _ ^^^^ (inhomogeneous), (14) 

in any dimension of space (including the dimension d = 1). This linear scaling has been observed numerically in (STj 
for spatially inhomogeneous ID self-gravitating systems, whereas for spatially homogeneous ID plasmas [H, [s^ the 
relaxation time scales as N'^to as explained previously. 

In certain cases, one can implement a local approximation and compute the collision term as if the system were 
spatially homogeneous. In that case, one obtains the Vlasov-Lcnard-Balcscu equation 

where <I'(r,i) = / u{\y — r' \) p{Y\t) dr' is the self-consistent potential created by the particles. In this equation, the 
effects of spatial inhomogeneity are retained only in the advection term (l.h.s.), while the collision term (r.h.s.) is 
calculated as if the system were spatially homogeneous. If we neglect collisions, we get the Vlasov equation which is 
valid for N — )■ -|-oo. If we neglect collective effects, replacing |e(k, k • v)| by 1, we obtain the Vlasov-Landau equation 

which may also be written as 

df 9f„ df d f , .w'^SiT — WiW.; f d d \ ,^ , ^ 
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We note that, for d = 1, the colhsion term vanishes so that the kinetic equations based on the local approximation 
reduce to the Vlasov equation. This would imply a relaxation time > Nt^ . Since the rigorous treatment of spatial 
inhomogeneity leads to a relaxation time ^ Ntjj (see above), we conclude that the local approximation is not valid 
in c? = 1 (see Paper III). 

Remark: If a ID system is initially spatially homogeneous but, evolving under finite N effects (collisions), becomes 
Vlasov unstable and undergoes a dynamical phase transition towards a spatially inhomogeneous distribution, we 
expect that the relaxation time will scale like 

^whoie _ js^^to, {d=l dynamical phase transition), (18) 

with 1 < 5 < 2 intermediate between Eqs. (fT3|) and (|14p . Such an anomalous scaling with 5 = 1.7 has been 
observed for the HMF model [ssj in the situation where the system undergoes a dynamical phase transition from a 
non- magnetized state to a magnetized state [soj . 



IV. RELAXATION OF A TEST PARTICLE IN A BATH: THE FOKKER-PLANCK EQUATION 

We now consider the relaxation of a "test" particle (tagged particle) in a bath of "field" particles with a steady dis- 
tribution function /(v). We assume that the field particles are either (i) at statistical equilibrium with the Boltzmann 
distribution (thermal bath), in which case their distribution does not change at all with time, or (ii) in the case of 
one-dimensional systems, in any stable steady state of the Vlasov equation (as we have just seen, in ci = 1, this profile 
does not change on a timescale of order Nto). We study how the test particle progressively acquires the distribution 
of the bath due to "collisions" with the field particles. The test particle has a stochastic motion. The evolution of the 
probability density P(v, t) that the test particle has a velocity v at time t is governed by a Fokker-Planck equation 
involving a diffusion term and a friction term. 



A. General expressions of the diffusion and friction coefficients 



The evolution of -P(v, t) can be obtained from the Lenard-Balescu equation ([3]) by considering that the distribution 
of the field particles is fixed (see Paper I for more details). Thus, we replace the time dependent distribution f{-v',t) 
by the static distribution /(v') of the bath. This procedure transforms the integro-differcntial equation ([3]) into the 
differential equation^'^ 



dk dv' ki kj 



u{kf 



e(k,k.v) 



r<5[k . (v - v')] 



_9_ 



d 



^(v,^)/(v'), 



(19) 



where e(k, uj) is the dielectric function corresponding to the fixed distribution function /(v). In fact, we can understand 
this result in the following manner. Equations ([3]) and (jl9p govern the evolution of the distribution function of a test 
particle (described by the coordinate v) interacting with field particles (described by the running coordinate v'). In 
Eq. (|3l), all the particles are equivalent so the distribution of the field particles /(v', t) changes with time exactly like 
the distribution of the test particle /(v, t). In Eq. the test particle and the field particles are not equivalent since 
the field particles form a "bath". The field particles have a steady (given) distribution /(v') while the distribution 
of the test particle /(v,i) = NmP{v,t) changes with time. In the BBGKY hierarchy, this amounts to singling out a 
test particle and assuming that the distribution function of the other particles is fixed. 
Equation ([T9| can be written in the form of a Fokker-Planck equation 



dP 
'dt 



d_ 

dv. 



dP 



-PF- 



pol 



involving a diffusion tensor 



= 7r(27r) 



dk dv' ki kj 



u{kf 
e(k,k.v)|' 



r<5[k.(v-v')]/(v'). 



(20) 



(21) 



The same result can be obtained by starting from the Fokker-Planck equation II25I I and calculating the first and second moments of the 
increments of velocity (see Paper I). In plasma physics, this Fokker-Planck approach, taking collective effects into account, was first 
performed by Hubbard [gsII . 
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and a friction force due to the polarization 

^2 



Fr'^A2nrm J dkdv' fc.fc, ^^j^^^<5[k • (v - v')]^(v'). (22) 



Using Eq. (PT|) . we can easily establish the identity 

dD 



dvj dvi ' 

where D — Da. If the velocity distribution of the field particles is isotropic, the diffusion tensor may be written as 



(23) 



D,,{v) = 



-y + j-jD^iv)S,,, (24) 



where I?||(f) and D±{v) are the diffusion coefficients in the directions parallel and perpendicular to the velocity of 
the test particle. We note that D = Da = + D±. 

Since the diffusion tensor depends on the velocity v of the test particle, it is useful to rewrite Eq. ([20)) in a form 
that is fully consistent with the general Fokker-Planck equation 



"J 
with 



^ = ^ (A P) - - 
dt dvidvn ^ dvi 



(AjP) - _^(Pi^/"=*-"), (25) 



^'^= 2At ' (^^^ 



The total friction force F friction is related to the friction due to the polarization Fpoi by 

Tpfriction j^pol , dDij 

Substituting Eqs. ([?!]) and ([^^ in Eq. (P7|) . and using an integration by parts, we find that the total friction force 
may be written as 

F/"''""' = ^(2^)^m / dkdV hk.fiy') ( ^ - ±-\ 5[k • (v - v')]^#L. (28) 



dvj dv'j J |e(k, k • v)P 



The two expressions ([20)) and ([25)) have their own interest. The expression ([25)) where the diffusion tensor is placed 
after the second derivative (DP) involves the total (dynamical) friction F friction and the expression pop where the 
diffusion tensor is placed between the derivatives dDdP isolates the part of the friction Ypoi due to the polarization. As 
shown in (90l - [9^ , the friction by polarization FpoZ arises from the retroaction of the field particles to the perturbation 
caused by the test particle. It can be obtained from a linear response theory. It represents, however, only one 
component of the dynamical friction ¥ friction, the other component being djDij (see Paper I). 
If we neglect collective effects, the diffusion tensor and the friction force reduce to 

= n{27rfm J dkdV k,kju{kf5[\i ■ (v - v')]/(v'), (29) 

(30) 



Fi friction 



Ff°^ = 7r(27r)'*77i j dkdV k,k,u{k)^S[k ■ (v - v')]^(v'), 
7:{27rfm J dkdV hk,f{w') {J- - <5[k • (v - ^r')]u{kf . (31) 



After a series of elementary calculations, we find that 

dD^ = (32) 

Combining Eqs. ([27)) and ([32)) . we obtain 



friction — poh 



2Fpoi, (33) 



so that the friction force F friction is twice the friction due to the polarization Fp^i (if the test particle has a mass m 
and the field particles a mass ruf, the factor 2 is replaced by the ratio (m + mf)/m (see Appendix [C)) . 
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B. Thermal bath: Maxwell-Boltzmann distribution 



We consider a thermal bath that is formed by field particles at statistical equilibrium having the isothermal 
(Maxwell-Boltzmann) distribution 

/(v)=(^)''^e-i''™^^ (34) 

Substituting the identity df /dv ~ —/3mf{v)v in Eq. (p2|) . using the (5-function to replace k-v' by k-v, and comparing 
the resulting expression with Eqs. (|2ip and (IMl) . we find that 

F[°' = -Dij{v)Pmvj = -D\\ {v)Pmv^. (35) 

The friction coefiicient ^ = D^^I3m is given by an Einstein relation expressing the fluctuation-dissipation theorem. We 
stress that the Einstein relation is valid for the friction due to the polarization Fpoi, not for the total (dynamical) 
friction F friction that has a more complex expression due to the term djDij. We do not have this subtlety for the usual 
Brownian motion for which the diffusion coefficient is constant. For a thermal bath, using Eq. ([55]) . the Fokker-Planck 
equation ([^0]) takes the form 



dP_ d_ 

dt dvi 



dP 

A,(v)( — +/3mPi;, 



(36) 



where DjAy) is given by Eq. (PT|) with Eq. p4p . This equation is similar to the Kramers equation in Brownian 
theory [93|. However, in the present case, the diffusion coefficient is anisotropic and depends on the velocity of the 
test particle. For t — >■ -|-oo, the distribution of the test particle relaxes towards the Maxwell-Boltzmann distribution 
Pg(v) = (/3m/27r)''/2e-''""'/2 of the bath. Since the Fokker-Planck equation (pS)) is vahd at the order 1/iV, the 
relaxation time scales as i^*'' ^ Ntn in any dimension of space (sec Sec. lIVPp . If the distribution P{\-,t) is 
spherically symmetric, using DijVj = D^\(v)vi, we can write the Fokker-Planck equation as 



dP _ 1 d 
dt v'^^^ dv 



dv 

For the Maxwell-Boltzmann distribution the diffusion tensor (PT|) is given by (see Paper I): 



(37) 



A.(v) = .i2.r ('^V'^- / dk'-^TTT^^e-^-'^. (38) 



27T J J k |e(k,k- v)|2 

Separating the integration over the modulus of k from the integration over its orientation, and setting k = k/fc, we 
obtain (see Paper I): 



where 



with 



Aj(v) = / d]^hkje-^P"'^^-^^ 2?^ W k . V , (39) 



V4x) = "^-^^ dfc, (40) 

2(27r)'^-5 n Jq [1 ~ ij{k)B{x)f + riikyC{x)^ 



r]{k) = -{277 fu{k)pmp, B{x) = 1-2x6-'=^ e^^dy, C(a;) = v^se"^'. (41) 



In these expressions, n = p/m \s the numerical density of particles and Vm = {l/PmY/'^ is the root mean square 
(r.m.s.) velocity in one direction. 
For ID systems, the Fokker-Planck equation reduces to 
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with a diffusion coefficient 



(43) 



where 



1 w-^ 



3 /•+00 



r]{k)^k 



2(27r)5 n Jo [1 - r]{k)B{x)Y + r]{k)^C{x)^ 



■dk. 



(44) 



The friction by polarization satisfies the Einstein relation Fpoi{v) = ^D{v)Pmv. 

When collective effects are neglected, the diffusion coefficient is obtained from Eqs. ([M)) and (|40p by taking 
B = C = 0. In that case, we obtain 



(45) 



where 



1 



3 /-H-oo 



2(27r)'^-5 n 7o 



Explicit expressions of the diffusion tensor can be obtained in various dimensions of space 0, H^l ■ 
In d = 3, we have 



(46) 



Aj(v) 



ViVi 1 , , „ 



(47) 



with 



G 



Dn{v) = 27r3/2l?3 



erf 



v] - G 



(48) 



G{w) 



2 1 



2w^ 



eri(w) r=e 



where erf(ii;) = JJ" e " is the error function. We note the asymptotic behaviors 



£'±(0) = ^—, i:)||(w)~+oo 



They imply Aj (v) ~ Z?||(0)% = ^D^{0)6,j when |v| ^ 0. 
In d = 2, we have 



D±{v) 



+00 



27r3/2p3 



A,(v) = [i?||(«) - D^{v)] ^ + i5^(z;)J.„ 



with 



h[—v 



D±{v) — 2?27re 



r , 2 



r ('^"^ 2 
h{—v 



where /„(cc) are the modified Bessel functions. We note the asymptotic behaviors 



Z?|| (0) = ttV2 , (0) = 7rl?2 , ^11 (v) 



, Dj_{v) 



20FP2 



(49) 



(50) 



(51) 



(52) 



(53) 
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They imply A,(v)~Di[(0)% 
In d = 1, we have 



D±{Q)S^J when |v| ^ 0. 



D{v) ^ 2Vie 



Using these analytical results, we can expUcitly check that formula ([5^ is satisfied. 

In the Dcbye-Hiickel approximation, the previous expressions remain vahd provided that is replaced by 



1 



-\-oo 



riDHikyk'^ dk, with r]DH{k) ^ -{2T:YuDH{k)l3mp 



yjk) 

l^vik)' 



(54) 



(55) 



This amounts to taking B ^1 and C = in Eqs. ([39]) and (|40)) . 

For spatially inhomogeneous systems, the Fokker-Planck equation must be written with angle-action variables as 
in |76l - l79l |. When a local approximation can be implemented, the distribution P{r, v, t) of a test particle in a thermal 
bath is governed by the Fokker-Planck equation 



dP dP , dP d 

h V • V<i>(r) • = 

dt dr dv dvi 



dP 

Aj(r,v) ( — +/3mPvj 



(56) 



where $(r) = J u{\r — r'|)p(r') dr' is the static potential produced by the field particles at statistical equilibrium and 
Dij{r,v) is given by Eqs. (pn|) - (|41l) where n is replaced by n{r). For ID systems, the foregoing equation reduces to 



dP 
'dt 



dP 
dx 



dP 
dv 



d_ 

dv 



D{x,v) 



dP 
dv 



PmPv I 



(57) 



C. Out-of-equilibrium bath in d = 1 

The derivation of the Fokker-Planck equation (jl9p . relying on a bath approximation, assumes that the distribution 
of the field particles is "frozen" so that /(v) does not change in time, at least on the timescale Ntn corresponding to 
the relaxation time of a test particle in a bath. This is always true for a thermal bath, corresponding to a distribution 
at statistical equilibrium (Boltzmann), because it does not change at all. For d > 1, due to collisions, any other 
distribution function relaxes towards the Boltzmann distribution on a timescale Ntu. Therefore, for d > 1, only the 
Boltzmann distribution can form a bath (the other distributions change in time). The situation is different in d = 1 
since the Lenard-Balescu collision term cancels out so that a Vlasov stable distribution function does not change on 
a timescale of order Ntu (see Sec. IIIIBp . Therefore, any Vlasov stable distribution function can be considered as 
"frozen" on this timescale and forms a bath. We can therefore develop a test particle approach in an out-of-equilibrium 
bath characterized by an arbitrary Vlasov stable distribution /(w), which is not necessarily the Maxwell-Boltzmann 
distribution of statistical equilibrium. 

For ID systems, the expressions (HU, and ([27)) of the diffusion coefficient and friction force take the form 



D{v) = 4Tr'^mf{v) 



-\-oo 



kujkf 
E(fc, kv)\- 



■ dk, 



(58) 



Fj,oi{v)^ATi''mf{v) 



+ CX3 



ku{kf 
\e{k,kv)\- 



■ dk. 



FfricUon{v) = FpolXv) + D' (v) , 



(59) 
(60) 



with (see Paper I): 



|e(fc,MI' 



f'(u) 

1 - 2nu{k)V I ^^-^ du 



+ 47r^u{kff{vf 



(61) 



We note the relation 



Fpoiiv)^D{v) 



d\nf 
dv ' 



(62) 
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which can be considered as a generalized Einstein relation valid for an out-of-equilibrium bath. The Fokker-Planck 
equation (jl9p may be written in the form 



dP_d_ 
dt dv 



fdP din/ 



(63) 



Equation (|63p is similar to a Fokker-Planck equation describing the motion of a Brownian particle in a potential 
U{v) — —\nf{v) created by the field particles. The distribution function of the test particle P{v,t) relaxes towards 
the distribution of the bath f{v) for t +00. Since the Fokker-Planck equation (|63l) is valid at the order 1/A^, the 
relaxation time scales as t^*^ ^ Nto- The fact that the distribution P{v,t) of a test particle relaxes towards any 
stable distribution f{v) of the field particles in d = 1 explains why the distribution of the field particles does not 
change on a timescale ^ Nto- As we have said, the situation is different in d > 1 [96j . 

For an isothermal distribution (thermal bath), we recover the results of Sec. IIVBI For the waterbag distribution, 
using Eq. (B.3) of Paper I, we obtain 







-| . 27ru{k)p 



(64) 



for —Vm V < Vm and D{v) — otherwise. 

When collective effects are neglected, the diffusion coefficient and the friction by polarization reduce to 



+00 

2 



D{v)^mf{v)l A-K^uikfkdk, (65) 





Fpoiiv) =mf'{v) 4:TT^u{kfkdk. (66) 

Jo 

They are respectively proportional to the distribution function f{v) and to its derivative f'{v). We immediately check 
that D'{v) = Fpoi{v) and Ffriction{v) = '2Fpoi{v) in agreement with Eqs. (|32l) and (pSj) . 

D. The relaxation time 

The relaxation time of a test particle in a bath may be estimated by i^*^ ~ ~ 1/(D/3to) ~ 1/C where D is 

the typical value of the diffusion coefficient. Using Eq. (|46p to estimate D, we obtain 



^bath or, equivalently, d"*'' n^;^ • (67) 

v^jt v{kYkddk nm^ ui^kfkd dk ^ ' 



We can take collective effects into account in the Debye-Hiickel approximation by simply replacing uik) by iiDHik), 
or r]{k) by rjonik). Using 77 ~ 1, ~ l/i?, n ^ N / R"^ and to ^ R/vm, we get 



i^j"*" - NtD, (Vd). (68) 

This is the timescale controlling the relaxation of a test particle in a bath, i.e. the time needed by a test particle to 
acquire the distribution of the bath. The timescale t^*'' should not be confused with the timescale t™''°''= controlling 
the relaxation of the system as a whole (see Sec. In d > 1 they are equivalent {t'^°^'^ ~ ^/j*'') but for spatially 
homogeneous systems in d = 1 they are totally different (^^''"'"^ > t^*'') since > Nto- In d = 1 dimension, the 

distribution P(w, t) of a test particle relaxes towards the distribution of the bath f{v) (that can be different from the 
Maxwellian as long as it is a stable steady state of the Vlasov equation) on a timescale of order Nto while the overall 
distribution of the system f{v, t) does not change at all on this timescale. For one dimensional plasmas, the difference 
in behavior of "distinguished" particles that relax on a timescale NtD and the overall population that relaxes on a 
timescale N'^tjj is shown in [sil. |86|. As a clear sign of the inequivalence of these two descriptions (evolution of the 
system as a whole and relaxation of a test particle in a bath), we note that the Lenard-Balescu equation ^ conserves 
the energy while the Fokker-Planck equation does not. 
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V. A SIMPLE APPROXIMATE KINETIC EQUATION WHICH DOES NOT DISPLAY ANY 
DIVERGENCE FOR 3D COULOMBIAN PLASMAS 

Prior to Lenard and Balescu, several authors [ssl - IH^ have attempted to derive from the BBGKY hierarchy a kinetic 
equation which does not display any divergence in order to improve the Landau theory. However, their treatment is 
approximate because, at some point of their derivation, they replaced some correlation functions by their equilibrium 
values (obtained with the Boltzmann distribution) while Lenard and Balescu have been able to determine the exact 
out-of-equilibrium correlation functions. Nevertheless, the approaches of [58l - l62j have the advantage of the simplicity 
and they clearly show how the divergences in the Landau equation can be eliminated. In this Section, we synthesize and 
generalize these different works and derive a simple approximate kinetic equation that does not display any divergence 
for 3D Coulombian plasmas. In the dominant approximation, it gives the same result as the Lenard-Balescu equation. 
We show, however, that this approximate equation breaks down in lower dimensions. 

The first two equations of the BBGKY hierarchy are 

^(1) = -1A. f F(2^1)g(l,2)dr2dv2, (69) 
at m ovi J 



rfrsrfvg^- 1^(1)5(2,3) + (lo 2), (70) 

OTi OVi 

where /(I) = /(vi,t) is the distribution function and g(l,2) — g(|ri — r2|,vi,V2,t) is the two-body correlation 
function. On the other hand, F(2 1) = ^ "^^j^ is the force by unit of mass created by particle 2 on particle 
1. These equations arc valid at the order 1/N. In particular, we have neglected the three-body correlation function 
which is of order l/N'^. This allows us to close the hierarchy. We note that the third term in Eq. ([70)1 is of order 
1/iV^ (see the scalings of Sec. |ll|) but this term becomes important at small scales when |r2 — ri| — > 0. Indeed, it 
is the term that accounts for strong collisions. We refer to 
notations as here^^ and for additional discussion. 

At statistical equilibrium, we expect the two-body correlation function to be of the form 



30| for a derivation of these equations using the same 



.g(|ri - r2|,vi, V2) = ./m(vi)/ji/(v2 )■(/;( In - r2|), (71) 
where /m is the Maxwellian distribution ([M)) . Substituting this ansatz in Eq. ([70)) . we get 

1^(1, 2) • VI + Pm'^ ■ vi^(l, 2) + (1^2) = -/3m2^ • vi - /3pm f dr,^ ■ viV(2, 3) + (1 o 2). (72) 

OTi OTi ori J OTi 

This equation being true for any vi and V2, we finally obtain 



^^-(l,2) + /3m^— ^(l,2) = -/3m^— -/3pm j dr3^V(2,3), (73) 

where x = ri — r2. This equation may also be derived from the YBG hierarchy starting from the microcanonical or 
canonical distributions [qt}. 

If we neglect in Eq. (|73|) strong collisions (second term) and collective effects (last term), we get after integration 

V'(x) = -I3m^u{yi). (74) 

In the case of plasmas, this expression is valid for Ai, ^ |x| <C \d where Al is the Landau length and A^) is the Debye 
length (see Paper III). This equation suggests interpreting Ue// ~ —%p/{f3m?) as an effective potential. 
If we only neglect strong collisions, Eq. ([75]) becomes 

J^(l,2) = _/3m2|!^-/3pm /dr3^^(2,3). (75) 



9x d. 



'X 



Actually, the equations in [30| are more general since they take into account effects of spatial inhomogcneity. They reduce to Eqs. I|69| l 
and I I70I I for spatially homogeneous systems. 
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Since the integral is a product of convolution, this equation can be easily solved in Fourier space leading to 



In the case of plasmas, this expression is valid for |x| ^ X^. This is the proper generalization of the Debye-Hiickel 
theory (the effective potential Ue// = —tp/lfSni^) is precisely the Debye potential udh given by Eq. (jA4p ). 
If we only neglect collective effects, Eq. ([73)) becomes 



^(1, 2) + /3m2|^^(l, 2) = ^Pm'^. (77) 
ox ox ox 

This equation can be easily integrated in physical space leading to 

i/;(x) = e-'^™'"^^) - 1. (78) 

In the case of plasmas, this expression is valid for |x| <C \d- Using Eqs. ([7T|) and ([75]) . we find that, when |ri — — > 0, 
the two-body distribution function of the A^-body problem defined by N^m'^P2{l,2) = /j\/(l)/j\/(2) [1 + x(li2)] is 
given by 

P2(ri,Vi,r2,V2) = ^-L_e-^^^("-------=), (79) 

where H2 = ^mvf + ^mvf + m^M(|ri — r2|) is the Hamiltonian of the two-body problem. Therefore, at small scales, 
the two-body distribution function is given by the canonical distribution of a pair of particles in interaction^^. For 3D 
gravitational systems for which u = —G/r, we find that V'(x) = —1 + e^^™" /'"L This correlation function accounts 
for the formation of "binaries" when the stars come at a distance of the order of the Landau length = G(5m? 
(the stars attract each other and form a pair). For one-component 3D plasmas with a neutralizing background for 
which u = / {m^r), we find that il>{yi) = — 1 + e^^'^"/!''' . This correlation function accounts for the repulsion between 
two like-sign charges when they come at a distance of the order of the Landau length Al = /3e^ (like-sign charges 
repulse each other). If we consider a more realistic two-components plasma made of positive and negative charges, 

the correlation between positive and negative charges is given by (x) = — 1 + e*^*^ ^^^K This correlation function 

accounts for the formation of "atoms" when opposite charges come at a distance of the order of the Landau length 
Xl ~ /3e^ (opposite-sign charges attract each other and form a pair). 
For 3D plasmas, Eqs. (1741) . ([76| and ([78]) may be written as 



^(x) = -^, ^(x) = -^e-l^l/^°, ^(x) = -l + e-^^'/l''l. (80) 
1^1 1^1 

These equations, that are valid at different scales, may be unified in a single equation 



V'(x) = exp 



-e 



1. (81) 



This expression was derived in |10d . llOlj in a different manner. 

Following [ssI - Ib^ , we can use these equilibrium results to solve Eq. (|70p approximately. The idea is to replace the 
correlation function g{\ri — r2|, vi, V2, t) in the third (strong collisions) and fifth (collective effects) terms of Eq. ([70]) 
by the approximate expression 

5(|ri - r2|,Vi,V2,i) ~ /(vi, i).f (v2, i)VX|ri - r2|), (82) 

where f(y,t) is the out-of-equilibrium distribution function and V'(|ri — r2|) is the equilibrium correlation function 
satisfying Eq. (j73p . After simplification, we obtain 

i§^(l, 2) + VI . 1^(1, 2) + (1 o 2) = • ^/(l)/(2) + (1 O 2), (83) 

2 ot OTi ori ovi 



This result is often assumed in the hterature (see, e.g., [9^ l99l|'l. ft is interesting to note that it can bo derived from the BBGKY 
hierarchy or from the YBG hierarchy. 
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where 

'0(X 



is an effective potential determined by the equihbrium two-body correlation function. Eqs. (|69p and (j83|) are now 
identical to those leading to the Landau equation [see, e.g., Eqs. (18) and (19) in [s^l] except that u is replaced by 
Me// in Eq. We can therefore directly write down the kinetic equation satisfied by /(v,t) as 



% = ^(2^)'™^ / dkdV hk,uik)u,jf{mk ■ (v - v')] - ^) /(v, t)/(v'. 



t). (85) 



This equation can also be written as Eq. © with replaced by K^^-^ defined with u{k)ueff{k) in place of u{kY . If 
we neglect strong collisions, we can use the Debye-Hiickcl expression (|76p of the correlation function, i.e. Ueff{k) ~ 
UDnik), and we obtain 

We stress that Eq. ([55]) is different from the Debye-Hiickcl approximation of the Lenard-Balescu equation [see Eq. 
([8])] since it involves u{k)uDH{k) instead of U£iH{kY ■ Actually, if we replace in the Lenard-Balescu equation ([3]) the 
quantity l/|e(k, k- v)]^by its value averaged over the velocity with the Maxwell distribution (f34|) . and use the identity 
(see Appendix C of fl02l|): 

/(V) ... ' ' ,8„ 



[2/ pj |e(k,k-v)|2 e(k,0) 1 + l3pm{2'KYu{ky 

we get the approximate kinetic equation (|86p . This gives a clear interpretation to that equation. 

In the thermal bath approach of Sec. IIVBI Eq. ([55)1 may be converted into the Fokker-Planck equation with 
a diffusion tensor given by Eq. (|45p with 

~ 2(2^)'i-V2 „, l-r;(fc) 
For a Coulombian plasma, using ri{k) = —kfy/k'^ (sec Paper III), we get 

25^^' = / dk. (89) 



2(27r)rf-i/2 n fc^ + /,2^ 

In d = 3, the integral in Eq. is convergent for fc — 0. Therefore, the approximate treatment of this Section is 
sufficient to suppress the divergence at large scales that appears in the Landau equation. On the other hand, the 
integral in Eq. (|89)) displays a divergence at small scales because we have neglected strong collisions in Eq. (|86p . We 
can cure this divergence heuristically by introducing a cut-off at the Landau length^^ . Had we taken into account the 
complete expression of the effective potential u^f / (r) this cut-off would not have been necessary since the correlation 
function (j78p accounts for a strong repulsion of like-sign charges when they come at a distance of the order of the 
Landau length^^. Indeed, the effective potential behaves like Ueff{k) ~ u{k)/{k/kL) for k — > -|-oo and the small scale 
divergence is suppressed. When strong collisions are taken into account, using 

u(k) 

Ueff{k) = ^f^, (91) 

1 + 1? + ^ 



^® The value of the resulting integral is 



L 



k-2 + kl 2 



dfc = i In (1 + A2) , (90) 



where A = k^/ko- 

^"^ For two-components plasmas and for self-gravitating systems, instead of a repulsion, we have an attraction between opposite-sign charges 
or between stars. This leads to the formation of "atoms" (-|-, — ) or "binary stars" which requires a special treatment. 
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Eq. (|89|) may be replaced by 



3 2(27r)5/2 n ""J^ k^ + kl + ^ ^ ' 



3 /■+00 



In the dominant approximation, we obtain 



In A 

2(2^)5/2 -m-^ — 



= 7:77^<kD-^, (93) 



where A = A_d/Al = is the number of charges in the Debye sphere. Similarly, K^^^ ~ (27re'*/m'^) In A. This 
is the same result as the one obtained from the Lenard-Balescu equation in the dominant approximation (see, e.g., 
Paper III). In d = 2 the integral in Eq. (|89p does not present any divergence (implying that strong collisions are 
negligible) and we get 

^^^'■' = 4(2^^™^-!' 

where A = nk^ . Similarly, K2^^ ~ ire^ /m^kjj. This slightly differs from the exact result obtained from the Lenard- 
Balescu equation (see Paper III). However, the approximate treatment of this Section suppresses the linear divergence 
at large scales present in the 2D Landau equation. Finally, in d = 1, the integral in Eq. ([M)) displays a logarithmic 
divergence at large scales that is absent in the Lenard-Balescu equation (see Paper III)^^. However, this logarithmic 
divergence is less severe than the quadratic divergence at large scales present in the ID Landau equation. 

As a final remark, we note that we can regularize the divergence at small scales in the Lenard-Balescu equation for 
3D plasmas by writing the kinetic equation ^ in the form 

I - '(^'f'r'i-J if^* ■ (V - v')l (A - ^) /(v,0/(v',t), (95) 

where x(fc) — > 1 for fc and x{k) for k +oo. The function x{k), which is related to il^ik), takes into account 
the strong correlations between the particles that appear at small scales. If we use the same arguments as those leading 
to Eq. (|?T|). we get xi^) ~ k^/k for fc — > -|-oo. In that case, Eq. (^5]) is convergent at both small and large scales. 
In the dominant approximation, we recover Eq. (j93p . If we assume that the particles are hard spheres^^, we may 
take xik) ~ l/(afc) for k — -|-oo where a is the size the particles. Assuming a small, and considering the dominant 
approximation, we get Eq. (|93p where In A is replaced by ln(A£)/a). More generally, we should have ln(A£)/A„im) 
where Xmin = max{AL,a} (in the gravitational case, we obtain the same results with the Debye length replaced by 
the Jeans length or by the size of the system). 



VI. THE EFFECT OF AN EXTERNAL STOCHASTIC FORCING 



The effect of an exter nal stocha stic forcing on the evolution of a system with long-range interactions has been 
studied by Nardini et al. |l03l Il04j and, independently, by ourselves in Paper I. Our approaches are different and our 
results also differ at first sight. In this section, we would like to discuss further this interesting problem and point out 
some connection between the two works. 

Let us first consider the case of particles without interaction submitted to an external stochastic force Fe(r, t). The 
equations of motion are 

^=v„ ^^-ev, + Fe(r,(t),t), (96) 

where we have assumed that the particles experience a friction force — with an inert medium. These equations are 
similar to the Langevin equations describing the Brownian motion of colloidal particles. However, in the standard 



This divergence only occurs in the thermal bath approach. In d = 1, the kinetic equation l|85|l . like the Lenard-Balescu equation, reduces 
to df/dt = 0. 

In general, the size of the particles is completely negligible in plasmas and stellar systems. A situation where the size of the particles 
matters is the formation of planetesimals and planets by gravitational instability. 
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approach [931, the stochastic force is written as \/2D'Ri(t) where Ri(i) is a Gaussian white noise such that (Ri(i)) = 
and {R'^{t)Rj{t')) = 5ij5^,y5{t — t'). Therefore, the present problem is a httle different. As in Paper I, we introduce 
the discrete distribution function /^(r, v,t) = m(5(r — Ti{t))S{w — Vi{t)) and write the continuity equation as 

dfd , dfd dfd d \ /^■7^ 

Proceeding as in Paper I, we decompose the exact distribution function as fd = .f + Sf where / = (fd) is the smooth 
distribution function and Sf represents the fluctuations around it. For a spatially homogeneous system / = /(v,<). 
Substituting this decomposition in Eq. dHZ]), we get 

df dSf dSf ^ df ^ dSf d , d , 

at ot or ov ov av av 

Taking the local average of this equation, we obtain the equation of evolution of the smooth distribution function 

f = -|;-(W> + |;-(C/v). (99) 

Subtracting Eq. (|99p from Eq. (|98p and neglecting the nonlinear terms (quasilinear approximation), we obtain the 
equation for the fluctuations 

.(eW--F..f (100) 
ot Or Ov Ov 

This equation may be solved formally as 

^/(r,v,t) = -^ G(t,t-r)Fe(r,t-r)|^(v,t-T)dT, (101) 

where G{t,t — r) is the Green function associated with the operator in the left hand side of Eq. pOOp . and we have 
assumed that (5/ = at t = 0. From Eq. (|10ip . we obtain 

f\Ftir,t)FHr,t~T))f{^,t~T)dT. (102) 



dvj JO 

If the auto-correlation function of the external force is given by 

{F,'{r,t)FJ{r,t')) = 2D6,,S{t - t'), (103) 

we get 

(<5/F,) = -D^. (104) 
Substituting this result in Eq. (|99)) . we obtain the equation 

^z,2i + yvV (105) 



dt dv \ dv 

In order to recover the Maxwell distribution at statistical equilibrium, the Einstein relation ^ = Dj3m must hold 
between the friction coefficient, the diffusion coefScient and the temperature. We note that Eq. pOSp coincides with 
the Fokker-Planck equation that is usually obtained from the standard Langevin equations with the Gaussian white 
noise V2Z?Ri(<). However, to obtain Eq. (jlOSp . we have not used the usual formalism of Fokker-Planck equations 
(Kramers-Moyal expansion, Ito formula...). In terms of the function P{k) defined in Paper I, the correlation function 
of the stochastic force may be written as 



{F^ir, t)F^ir', t')) = 2^5^j5{t - t') J fc2p(fc)e^k.(r-r') ^j^^ ^^gg^ 
Comparing this equation with Eq. (|103p . we obtain the following expression of the diffusion coefficient 

D = Try k^P{k)dk. (107) 
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This expression is valid when there is no interaction between the particles. 

In Paper I, we have tried to generalize this approach in the case where the particles interact via long-range forces. 
In that case, wc must add the force —V^d{^i{t),t) = — J Wu{\ri{t) — r'\)fd{r',v',t)dr'dv' produced by the other 
particles in the equations of motion (|96p . The continuity equation then becomes 

?^+v.f^-V<I>..fi+F..?/^^ = A.(^/,v). (108) 
ot ar av ov ov 

Within the quasilinear approximation, wc obtain the closed coupled equations 

| = |,-(^/V.<^>)-|;-(W) + |;-(e/v), (109) 

and 

where we have assumed that the friction is weak (^ <^ 1) so that it can be neglected in Eq. (|110p . From these 
equations, we have obtained in Paper I a nonlinear Fokker-Planck equation of the form 

'4-l(^m^+ifv,). (Ill, 



dt dvi \ ' dvj 



with a diffusion tensor 



20 



In the non-interacting case, we have e(k, k • v ) = 1 and we recover Eq. (|107p . 

Using a different approach, Nardini et al. [l03l Il04j arrive at a different result (we can make contact with their 
notations by setting Ck = 27rfc^P(fc)). We briefly try to find out a connection between the two approaches. First of all, 
we note that when Fe(r, t) is a stochastic force, the validity of the continuity equation (|108l) is not obvious. Indeed, 
if we use the Ito formalism in the spirit of (lOSl Il06| , we get an equation of the form 

— + V • — - • — + F . ^ = A . (d— + ifdw] (113) 

where D is given by Eq. (jl07p . Then, within the quasilinear approximation, Eq. (jlOOp is replaced by 

while Eq. (|110p remains unchanged in the limit of weak friction s £ <g; 1. In order to obtain Eq. (|114p . we have treated 
the term /^Fe in Eq. (|113p as a noise (in the same spirit as in |l05l . Il06j ) and, consequently, we have not written the 
term (SfFe) in Eq. (|114p . In this way, in the absence of interaction, Eq. (|114p reduces to Eq. (|105p as it should. 
On the other hand, the term (SfWS^) has been computed in Appendix C of Paper I. Substituting this result in Eq. 
(|114p , this leads to a nonlinear Fokker-Planck equation of the form 

with a diffusion tensor 



dt dvi \ '-^ dvj J 9v \ 9v 



This equation coincides with the result in jl03l Il04l |. However, if we evaluate the term {6fFe) with the method of 
Paper I, we find that it is non-zero and that the sum of the terms (SfF^) and —{SfW6^) simplifies to give the diffusion 
coefficient (fTT2)) . 



When considering the effect of the external stochastic force, the Laplace transforms in Sec. 6 of Paper I must in fact bo replaced by 
Fourier transforms. As a result, Eqs. (141) and (C.4) of Paper I must be multiplied by (27r)'^. 
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VII. CONCLUSION 



111 the introduction of this paper, we have proposed a short historic concerning the development of kinetic theory in 
plasma physics. We have mentioned important contributions that are not well-known. In particular, we have noted 
that, prior to the classical works of Lenard and Balescu, some authors derived an approximate kinetic equation that 
solves the divergences at small and large scales appearing in the Landau equation. In the main part of the paper, 
we have developed these works further and we have shown that this approximate kinetic equation may be obtained 
from the Lenard-Balescu equation by replacing the quantity l/|e(k, k • v)p by its average (l/|ep|) over the velocities 
computed with the Maxwell distribution. We have also recapitulated the main kinetic equations of plasma physics 
and we have extended them to arbitrary long-range potentials of interaction in different dimensions of space. We 
have discussed the scaling of the relaxation time with the number of particles and with the dimension of space. We 
have also briefly discussed the effect of a stochastic forcing on the evolution of systems with long-range interactions. 
In Paper III, we apply these results to specific systems with long-range interactions such as systems with power-law 
potentials, plasmas and stellar systems in arbitrary dimensions of space, and the HMF model with attractive ad 
repulsive interactions. 



Appendix A: The Debye-Hiickel theory for a general long-range potential of interaction 

In this Appendix wc provide the proper generalization of the Debye-Hiickel theory for an arbitrary long-range 
potential of interaction. We wish to determine the effective interaction between particles at statistical equilibrium 
taking collective effects into account. A single (test) particle at rp = produces a "naked" potential ^o{r) = mw(r). 
This potential modifies the distribution of the other (field) particles. The resulting change of density p{r) in turn 
produces an extra-potential which adds to the original one 5'o(r). The effective, or "dressed", potential created by 
the test particle is therefore the solution of the self-consistency equation 

^effir) = J u{r - r') [p{r') + m5{r')] dr' . (Al) 

The variation of the distribution of the field particles is given by the Boltzmann statistics 

p{v) = pe-'3"*<=//('-) - p, (A2) 

where p is the uniform distribution of the particles. In the weak coupling limit ^m^eff ^ 1, we can expand the 
exponential in Eq. (jA2[) so that p ~ —j3mp'^(.ff- Substituting this result in Eq. (|Aip we obtain an integro-differential 
equation 

$e//(r) = I-') [-Pmp^effir') + m5{v')] dr' , (A3) 

determining the effective potential. Since the integral is a product of convolution, this equation can be solved easily 
in Fourier space. Writing $e/^(r) = mu£)fl'(r), we obtain 

UDH(k) = , , . (A4) 

Specializing this expression to 3D plasmas, we obtain the Debye-Hiickel potential 

or, in physical space, udh ~ (e^/TO^)r~^e~'^°'', where kr, = {'^■Ke^^plrfif-^'^ is the Debye wavenumber. An equivalent 
result, expressed in terms of the two-body correlation function, can be derived from the YBG hierarchy j97| or from 
the BBGKY hierarchy (see Sec. |V|. The generalization (IA4I) allows us to consider other potentials than just the 
Coulombian potential. Some explicit examples are given in |97| . 
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Appendix B: The Rosenbluth potentials in d dimensions 

In this Appendix, we show that the Landau equation may be expressed in terms of the Rosenbluth potentials [27| . 
Our results are valid in d dimensions. The Landau equation (jH]) may be written as 

'^-If '<„ i - /l^l .V, with . /O^ifc^, (Bl, 



dt dvi J \ dvj dii , ~ 
where / — /(v, t), f = /(v', t), and w = v — v'. It can be put in the Fokker-Planck form 



df ^ f ^■f fTpPol\ (^f /T^ f\ ^ f fj-, friction-. /T3r,\ 



where the diffusion and friction coefficients are given by 



3 ^ 

Using the identities 



= Kd „ „ , 
oviOVj 



and ^ = -(d- l)Kd^ = (d - l)Kd^ ^ ^ 
ovj w-^ ov, 



the coefficients of diffusion and friction may be rewritten as 



oviOv 



dh 



Ffr,cHon = 2Fpol = 2{d - l)/^d^(v,0, 



where 



5(v,t)= / /(v',t)|v-v'|dv', h{v,t) 



(B4) 



(B5) 



(B6) 



(B7) 



|v - v\ 

are the so-called Rosenbluth potentials [27[. In terms of these potentials, the Landau equation may be rewritten as 



^^Kd — 
dt dvi 



d^g df _ dh 
dvidvj dvj ' di'i 



df 
at 



dvidvj \dvidvj 



(B8) 



For an isotropic distribution function / ~ f{v,t), implying g ~ g{v,t) and h — h{v,t), the difl:usion tensor (|B6l -a) 
may be written as 



Dn 



d-1 



d-l 



with 



D\\ = K, 



D^ = {d-l)Kd--. 



On the other hand, the friction force ([B6]-b) may be written as 



1 dh 



F friction = 2Fpo; — 2{d — V. 



(B9) 



(BIO) 



(Bll) 



Using DijVj — the Landau equation (|B2[) becomes 
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Explicit expressions of g{v,t) and h{v,t) in terms of f{v,t) are given in [98[ in d = 3 (see also [30[). The Landau 
equation (|B12p describes the evolution of the system as a whole. This Landau-type equation, corresponding to a 
microcanonical description in which the energy is conserved, was used by King to study the dynamics of globular 
clusters. If we replace f{v,t) by P{v,t) and f{v',t) by the Maxwellian distribution f{v') given by Eq. ([M)) . we 
obtain the Fokker-Planck equation (|37p describing the relaxation of a test particle in a bath of field particles at 
statistical equilibrium (thermal bath). In that case, we recover the results of Sec. IIVBI This Fokker-Planck-type 
equation, corresponding to a canonical description in which the tem perature is fixed, was used by Chandrasekhar 
[isl [Toj to determine the evaporation rate of globular clusters (see [l07| for a more detailed discussion of the difference 
between the approaches of Chandrasekhar and King). Eq. (jB12p may be used to obtain a more general Fokker- 
Planck equation. Indeed, if we replace f{v,t) by P{v,t) and f{v',t) by any isotropic distribution f{v'), we obtain a 
Fokker-Planck equation describing the relaxation of a test particle in an out-of-equilibrium bath of field particles^^. 
We recall, however, that the results of this Appendix assume that the system is spatially homogeneous (or that a 
local approximation can be implemented) and neglects collective effects. 



Appendix C: Multi-species systems 

It is straightforward to generalize the kinetic theory for several species. The Lcnard-Balescu equation is replaced 

by 



with 



f k • ^ 

e(k, uj) = l + (2TTfu{k) y / ^ dv, (C2) 

^7 u; - k • V 

b 

where /°(v, t) is the distribution function of species a normalized such that J dw — Uama- We can use this equation 
to give a new interpretation of the test particle approach developed in Sec. IIVI We make three assumptions: (i) We 
assume that the system is composed of two types of particles, the test particles with mass m and the field particles 
with mass nif, (ii) we assume that the number of test particles is much lower than the number of field particles; 
(iii) we assume that the field particles are in a steady distribution /(v). Because of assumption (ii), the collisions 
between the field particles and the test particles do not alter the distribution of the field particles so that the field 
particles remain in their steady state. The collisions of the test particles among themselves are also negligible, so 
they only evolve due to collisions with the field particles. Therefore, if we call /'(v, t) the distribution function of the 
test particles (to have notations similar to those of Sec. IIVI with, however, a different interpretation), its evolution is 
given by the Fokker-Planck equation obtained from Eq. (|C1[) : 



with 



r k- ^ 

e(k, w) = 1 + {2TTfu{k) / ^ dv. (C4) 

J — k • V 

The diffusion tensor is given by 

A, - ^(2^)^m/ J dkdV fc,,A;,^^^-^-^5[k • (v - v')]/(v'), (C5) 



and the friction force due to the polarization by 



Fr'^n{2nrm J dkdV hk,j-^^^^5[k ■ - (C6) 



Actually, this more general approach is not self-consistent in d > 1 since, as already explained, the distribution f{v) is not steady on 
the timescale tn ~ Nt^, over which the test particle relaxes, unless f{v) is the Maxwell distribution. 
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We recall that the diffusion coefficient is due to the fluctuations of the force produced by the field particles, while the 
friction by polarization is due to the perturbation on the distribution of the field particles caused by the test particles. 
This explains the occurrence of the masses to/ and m in Eqs. (jC5[) and (jC6P respectively. 
If we neglect collective effects, using Eq. ([??)) and noting that 



V 



FF 



(C7) 



we get 



■ friction 



= 1 



TO/ 
TO 



pol • 



(C8) 



If we assume furthermore that m 3> to/, we find that F friction — Fpo;- However, in general, the friction force is 
different from the friction by polarization. The other results of Sec. IIVI can be easily generalized to multi-species 
systems. For multi-species systems, the dimension d = 1 is special as discussed in Sec. 7 of [s^ ]. 
In the thermal bath approach F[°'' = —Dij{\-)/3mvj and the Fokker-Planck equation reduces to 



dP 
'dt 



d_ 



dP 

A,(v)(— +/3TOPi;, 



(C9) 



where Dij{\-) is given by Eq. (|C5P determined by the field particles. If the distribution of the field particles is 
/ cx e~'^™^""/^. the equilibrium distribution of the test particle is P cx e-'3™"^/2 jm/nif ^ When to ^ to/, the 
evolution is dominated by frictional effects; when to ^ to/ it is dominated by diffusion. The relaxation time scales 
like t^*'* ~ v'^/D where v'^ is the r.m.s. velocity of the test particle (with mass to) and D is determined by the 
field particles (with mass to/). Therefore, f^*'* ^ l/{Dl3m) ^ l/S, where ^ is the friction coefficient associated to the 
friction by polarization (the "true" friction coefficient is = (1 + TO//m)^). 



Appendix D: Dimensionless Fokker-Planck equation 



With the change of variables w — (/3to/2)^/^v, the Fokker-Planck equation (p6|) may be rewritten as 



dP 1 d 



dt 2v^„ dwi 



A,(w) 



dP 



2Pw, 



(Dl) 



with 



A,(w) = / dkfc,fc,e^('"-)'2?d(k-w). 



(D2) 



When collective effects are neglected, using Dij{w) ~ 'DdGij{w) it may be written as 



dP _ 1 d 
dt tji dwi 



G,(w)(^ + 2P., 



(D3) 



with 



Gy(w) = / dkfeifcje-*-^)' 



tR 



2r? 



(D4) 



These scalings clearly show the emergence of the relaxation timescale f^. 
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